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The evolution of the dark matter distribution at the Galactic center is analyzed, which is caused by the 
combination of gravitational scattering on Galactic bulge stars and absorption by a supermassive black 
hole at the center of the bulge. Attention is focused on the boundary condition on the black hole. It is 
shown that its form depends on the energy of dark matter particles. The modified flux of dark matter 
particles onto the black hole is calculated. Estimates of the amount of dark matter absorbed show that 
the fraction of dark matter in the total mass of the black hole may be signihcant. The density of dark 
matter at the central part of the bulge is calculated. It is shown that recently observed 7 radiation 
from the Galactic center can be attributed to the annihilation of dark matter with this density. 
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1 Introduction 


Investigations of the distribution of dark matter in 
the nuclei of galaxies currently attract the interest of 
many researchers in view of both the effect of dark mat¬ 
ter on the growth of supermassive black holes and the 
search for its possible annihilation radiation. In this 
work, we analyze the evolution of the dark matter distri¬ 
bution at the Galactic center under the effects of scatter¬ 
ing by stars at the Galactic center and absorption by the 
central black hole. We assume that dark matter consists 
of particles that undergo only gravitational interaction 
and possibly slightly annihilate (so-called weakly inter¬ 
acting massive particles, WIMPs) p. 

The interaction of the central black hole with its en¬ 
vironment, including dark matter at the Galactic center 
(bulge), was analyzed in PlOiniQE!] and elsewhere. An 
approach in which the distribution function of the dark 
matter is written in terms of invariants of motion seems 
to be most consistent. In the case of spherically symmet¬ 
ric, sufficiently slow evolution (compared to dynamical 
time), such variables are the orbital angular momentum 
TO, its projection TOz onto the z axis, and the radial ac¬ 
tion (adiabatic invariant) 


I =- 

TT 


r^+ 


;(r) dr . 


We use a model of the dark matter halo formation in 
which its spatial density prohle before the formation 
of the Galaxy and baryonic compression had the form 
p{r) ~ and the initial distribution function in 

terms I, to, to^ has the form [3 

f{I,m,m„) = 5{m^-I q^I"^) , ( 1 ) 


where Zq ~ 0.1 is a small parameter. We emphasize the 
convenience of choosing radial action I rather than parti¬ 
cle energy E as a. basic variable when analyzing evolution 
in a slowly varying potential. 


1.1 Kinetic equation. In the zeroth approxima¬ 
tion, the distribution function of dark matter in the 
chosen variables does not change upon baryonic com¬ 
pression. Under the assumption that the black hole 
is formed at the coordinate origin, the increase in the 
black hole mass M^h causes the absorption of particles 
whose angular momenta to are less than the critical value 
TOg = AGMbh/c, but the amount of such particles turns 
out to be negligibly small 0. However, the situation is 
significantly changed due to effects associated with the 
gravitational scattering of dark matter particles on bulge 
stars. These effects lead to the diffusion of particles in 
phase space glEl. 

In the first approximation, the evolution of the dis¬ 
tribution function is described by the kinetic equation 


dfi{Ii},t) 

dt 


d 



( 2 ) 


where {R} = {/, to,to^}, and Rki are the corresponding 
diffusion coefficients. Owing to the spherical symmetry 
of the evolution, does not enter into the equations. 
In addition, I can be approximated as I{E, to) « J{E) — 
/3m [2], where /3 = I for the Goulomb potential of the 
black hole, and /3 ~ 0.6 for the isothermal potential of 
the bulge. Hereafter, we will use J instead of I, because 
we are interested in the region of low angular momenta. 
Since the parameter Iq is small, / can also be changed 
to J in the initial distribution O- It can be shown that 


I 





the cross terms R 12 are small in the region of interest 
13 . Finally, the kinetic equation takes the form 


df{J,m,t) 

dt 


A 



lA 

m dm 


m R 22 


A 

dm 


(3) 


The diffusion coefficients are calculated in [S]. One¬ 
dimensional diffusion along the m axis turns out to be 
the most significant. It leads to the particle flux into 
the region of low orbital angular momenta and to their 
absorption by the black hole for m < mg. This process 
was analyzed in detail in HIHl. The effect considered 
there is the following. Let us write Eq. m in the one¬ 
dimensional form 


df{J,m,t) d f df \ 

dt m dm \ ^ dm) 


(4) 


The initial condition is taken in the form CJ, and the 
boundary condition at the black hole corresponds to ab¬ 
sorption: /(J,TO = TOg,t) = 0. The diffusion coefficient 
R 22 calculated for the isothermal bulge has the form 


i ?22 = 0.46 (To, (5) 


1.2 Bulge model. We assume that the star dis¬ 
tribution in the bulge is isothermal, i.e., spherically sym¬ 
metric and isotropic, and has the power-law density pro¬ 
file 

( T \ —^ ^ 

/ ’ Po^o ~ 2 tt~G ’ 

where (Tq is the star velocity dispersion independent of 
the distance to the center. In the inner part of the bulge 
- the black hole influence region - the density profile has 
lower exponent [^: 

Ps{in){r) = Po , (10) 

where po and Rq are the same values as in ®, which 
assures the continuity of the star density at the boundary 
of the black hole influence region. We define the radius 
of the influence region such that the total mass of stars 
inside the region is equal to the mass of the black hole: 

Min = / Ps,in{r') iTTr'^ dr' = —poRo^^^ = Mbh , 

Jo 4 


where ag is the velocity dispersion of bulge stars. Mg 
is the mass of a star (for simplicity, we assume that all 
stars have the same mass equal to Mq), and Lc ~ 10 is 
the Coulomb logarithm. 

In this formulation, the solution of the diffusion equa¬ 
tion yields the following expression for the dark matter 
flux onto the black hole 


GMbh 4 2 


( 11 ) 


In Milky Way, (Tq = 85 90 km/s, and the value of Rg 

corresponding to the observed black hole mass M^h = 
3 • 10®M 0 , is Rg = 1.3 pc. These values also give the 
observation-consistent normalization for the star density 
profile at the center [See Ea. lUnil ]. 
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Here 
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is the flux in one-dimensional diffusion equation 0 with 
the initial condition = 0 ) = S{m'^ — (Igd)"^). 

However, the absorbing boundary condition is valid 
only for Ato^ <C where 


Ato = v' 2 r(J) R 22 ^ ( 8 ) 

is the rms change in the orbital angular momentum per 
orbital period T{J). It is easy to show that this condition 
is violated even at rather low energies of the particle. 
In this work, we obtain a more accurate expression for 
the boundary condition and investigate its effect on the 
absorption rate of the dark matter. We start with the 
presentation of the Milky Way bulge model. 


Thus, as was mentioned above, the absorbing-boun¬ 
dary approximation f{mg, J,t) = 0 is valid only if Ato ^ 
mg, i.e., change in the orbital angular momentum of the 
particle per period is small compared to the characteris¬ 
tic problem scale - the boundary orbital angular momen¬ 
tum. At the same time, as is easily seen from Eqs. (0 
01, this condition is violated for Milky Way at present, 
because the maximum of the flux comes from the values 
of J for which Ato ^ mg. 

To describe correctly the absorption of particles by 
the black hole, we consider two limiting cases (see fig¬ 
ure): the random-walk approximation (absorption for 
Ato ^ TOg, and the pinhole approximation for Am ^ 
TOg (the names are taken from Lightman&Shapiro | 1 ()|. 
where this process was considered in application to stars). 

The absorption models are different in the random- 
walk and pinhole approximations. In the former case, 
the distribution function f{m,J,t) 0 for to ^ TOg, 
and its derivative at the absorption boundary to = TOg 
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Figure 1: Two absorption regimes 


determines the flux. In the latter case, f(rn) —>■ fg where 
fg is a nonzero boundary limit, and the absorption rate 
is determined as follows. There is a nonzero probabil¬ 
ity per orbital period that the angular momentum m 
of a particle increases by 5m such that m -f 5m < mg 
and, hence, the particle is captured by the black hole. 
Flux onto the black hole is determined as the sum of 
the absorption probabilities for all particles per period. 
At the same time, it is equal to the diffusion flux from 
higher m values, where the diffusion approximation is 
valid (Am < m). 


2.1 Boundary condition. Different absorption 
regimes can be described by the following modification 
of the boundary condition: 


/ - mgU 
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dm 


= 0 . 


( 12 ) 


Values q : 1 correspond to the random-walk approxi¬ 

mation, and a 3> 1 to the pinhole approximation. 

Therefore, if we take the value of the distribution 
function at the boundary f(m = mg) = fg and assume 
that the flux is continuous near the boundary and varies 
only slightly at low m, then we obtain the following ex¬ 
pression for / near the boundary: 

1 777 

/(m)=/,(l-f-ln—) (13) 

Cx Iflg 

Our aim is to obtain expressions for a and for modified 
flux Sj{t) of particles with radial action J on the black 
hole in the limit a ^ 1. 

For this purpose, we calculate the number of parti¬ 
cles that diffuse into the region m < mg and are ab¬ 
sorbed by the black hole during one orbital period. We 
assume that the probability distribution for obtaining a 
given increment 5m of the particle angular momentum is 
Gaussian with the variance 5m^ = Am^ and mean value 


5m <C Am cni: 


p(m, 5m) 


1 / (5m^ \ 


1 -I- 5m jm 


l + erf(^^) 

'' v2 Am ' 


Am 


^exp( 2 An 


(14) 


The normalization is taken such that p(m, 5m) d5m = 

1 . 

The total number of particles absorbed per period is 
given by the expression 


poo p — {m—mg) 

Q= 2mdm/(m) / p{5m) d5m = Sj T{J). 

J rrifj J —m 

(15) 

At the same time, it is equal to the flux from larger m 
values per period: 


Q = SjT{J) =mR^T = ^RT = . (16) 

am a 2a 

In the limit Am ^ mg, we obtain 


^ ^ / Am 

a = 2.8 - 

V mg 


which agrees with the result obtained in nni. 


(17) 


2.2 Dark matter absorption. To determine the 
expression for the particle flux onto the black hole Sj{t), 
we use the same procedure as for the pinhole limit con¬ 
sidered in 1^. Namely, knowing the general form of the 
expression Sj{t) cx iexp(—|^), we combine the solu¬ 
tion /(m, t) given by Eq. I13II and the Gaussian diffusion 
solution for the initial condition f{t = 0) = 5{imf — mo^). 
Thus, we obtain 



This expression for flux differs from the one obtained 
in previous investigation o by an additional term a in 
the denominator. As was pointed out, at present a{J) ^ 
1 for the values of J from which flux is maximal and, 
therefore, the total flux is much lower than the value 
obtained without corrections. Indeed, the total flux is 
given by the expression 

poo 

S{t)=2{2n)^ dJ hJ^/^Sj{t)^ 

==9.5/o(^)''''“l-‘(/3+ir»'» (19) 
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where the quantity 


_ i/5 i ?22 1 27ri?22 ^ ^ f ^ ^ f ^ 

^ “ 20 /ocro 2 mg 2 VlOi° y ^3- lOeM©; 

( 20 ) 

reflects the correction due to the modihcation of bound¬ 
ary condition If the correction is disregarded (by 

setting /3 = 0 ), the expression flm yields a power-law 
increase of the black hole mass Mbh{t) cx [^. The 
estimate for the mass directly depends on the parameters 
/o and Iq of the initial dark matter distribution. Taking 
/o = 6-10® g (cmVs)-®/® and Iq = 0.1 for Milky Way ac¬ 
cording to in], we estimate the mass as Mbh — 11 • 10 ®Mq 
for the present time {t = 10^® yr = 3 • 10^"^ s) and for 
negligibly low seed black hole mass (under the assump¬ 
tion that mass increases only due to the absorption of 
dark matter; i.e., it is a lower estimate for the mass). 

The inclusion of the flux modihcation signihcantly 
changes the situation. As follows from Eq. 12011. if the 
black hole mass is low, then /3 ^ 1 and hence the dark 
matter hux is also low. This is an evidence of the bary- 
onic nature of a seed black hole. An upper estimate 
for the absorbed dark matter mass can be obtained 
as the difference = Mq — Mbh(O) between the H- 
nal mass Mbh{to = 3 • 10^”^s) = Mq = 3 • 10®Mq, and 
the initial mass MbhiO) obtained by solving the equa¬ 
tion dMbh/dt = S{t) backwards in time. For the values 
/o and Iq adopted above, we obtain Md = 0.67 Mq = 
2 • 1 O®M 0 . If the value /o is ten times lower, then 
Md = 0.11 Mq. Another estimate can be obtained by 
the supposition that the growth of the black hole due 
to absorption of both dark and baryonic matter has a 
power-law form Mbh{t) = Mo(f/fo)®'- In this case, the 
mass of dark matter absorbed is given by the expression 
Md = /q° S{t, Mbh{t)) dt. For the growth law with the 
exponent 7 = 1/2 (obtained in analyzing the absorp¬ 
tion of stars) Md = 0.36 A/q; it is proportional to the 
normalization constant /o. 

Thus, the correct formulation of the boundary condi¬ 
tion considerably (approximately by an order of magni¬ 
tude) reduces the estimate for the mass of dark matter 
absorbed by the black hole. The absorbed mass com¬ 
prises up to 1/3 — 1/2 of the current black hole mass for 
the chosen normalization of the dark matter density. We 
emphasize that a similar consideration for the absorp¬ 
tion of stars shows that the modification of the bound¬ 
ary condition is insignificant j^. This difference results 
from the fact that dark matter particles initially have 
low orbital angular momenta and sooner reach the ab¬ 
sorption boundary. This means that the flux at a given 
time comes from higher J values than those for stars. 
Hence, the boundary condition for dark matter changes 
its form at large J values more considerably. 


3 Density profile of dark matter and 

ITS ANNIHILATION 

We consider the detection of 7 radiation from the 
Galactic center, which is probably due to annihilation of 
weakly interacting particles of dark matter (WIMPs), as 
a possible observational test of the current dark matter 
distribution. Photons of energies of several TeVs have 
been detected on the H.E.S.S. telescope, whose angular 
resolution of 3' — 5.8' corresponds to a spatial region of 
about 10 pc in size (the distance to the Galactic center 
is taken to be 8.5 kpc). As was shown in im, the ob¬ 
served photon flux can be explained by the annihilation 
of supersymmetric particles with a mass of about 12 TeV 
whose density has a power-law profile and a mean value 
of about ~ IO^Mq/pc® inside the central 10 pc. 

Knowing the distribution function of dark matter in 
phase space, it is easy to calculate its density profile as 
a function of the distance from the coordinate origin: 


00 'f^max 


( 21 ) 

Let us calculate the dark matter density in the central 
region 10 pc in size. For simplicity, we assume that the 
potential in this region is determined by the central star 
cluster with a density Ps{r) = ao^j2TrGr^, where tro is 
the star velocity dispersion. As follows from P l ITH Il . 
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For r ^ 10 pc a ^ 100, therefore, the first fraction in 
Eq.J23) is close to unity. The quantity mmax in Eq. 
is determined from the condition of zero denominator: 
'^max = 2cro?'\/ln[0rJ/(Tor] ~ 2aQr for J > = 

o’o’'/v^- However, approximation 12211 is applicable only 
for low m values and, in particular, the total integral 
jrnmax cannot be larger than /o ac¬ 

cording to the normalization condition. This gives the 
restriction m < mo = V5 Rt, which becomes important 
for r > r* ~ 2 pc. Thus, the internal integral in Eq. en 
is given by the expression 
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r rrim ax = 2 cror , r < r* 

\ 2ao{r - - r* 2 ) , r > r* 

Since dE/dJ = 2 tt/T{J) = 2cto^/J, we have 
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In this work, the diffusion of dark matter at the 
Galactic center has been considered. This diffusion is 
caused by scattering of dark matter particles on bulge 
stars and is considered in the {J,m,mz} phase space, 
where J is the modified radial action. Diffusion along 
the axis of the orbital angular momentum to plays the 
main role in the absorption of dark matter. The pres¬ 
ence of the black hole determines the boundary condition 
at TO = rUg that has different forms for small and large 
values of J. The dark matter flux S{t) calculated using 
the modified boundary condition is much smaller than 
the value obtained in previous works, where this mod¬ 
ification was disregarded (i.e., the absorbing boundary 
approximation was applied). The amount of dark mat¬ 
ter absorbed by the black hole may at present compose a 
significant fraction of the black hole mass. The dark mat¬ 
ter density profile induced by diffusion within the central 
region 10 pc in size can be responsible for the observed 
7 radiation from the Galactic center, which arises upon 
the annihilation of dark matter particles. 
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